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Abstract 

We demonstrate the effectivity of the covariant background field 
method by means of an explicit calculation of the 3-loop /3-function for 
a pure Yang-Mills theory. To maintain manifest background invariance 
throughout our calculation, we stay in coordinate space and treat the 
background field non-perturbatively. In this way the presence of a back¬ 
ground field does not increase the number of vertices and leads to a rela¬ 
tively small number of vacuum graphs in the effective action. Restricting 
to a covariantly constant background field in Fock-Schwinger gauge per¬ 
mits explicit expansion of all quantum field propagators in powers of the 
field strength only. Hence, Feynman graphs are at most logarithmically 
divergent. At 2-loop order only a single Feynman graph without sub¬ 
divergences needs to be calculated. At 3-loop order 24 graphs remain. 
Insisting on manifest background gauge invariance at all stages of a cal¬ 
culation is thus shown to be a major labor saving device. All calculations 
were performed with Mathematica in view of its superior pattern matching 
capabilities. Finally, we describe briefly the extension of such covariant 
methods to the case of supergravity theories. 
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1 Introduction 


The essential quantity in quantum field theory is the effective action, which for 
a given classical action sums up all quantum corrections and from which the 
S-matrix can be found. Unfortunately, in practice this functional can only be 
determined in perturbation theory. In particular, for gauge theories the conven¬ 
tional effective action, in contradistinction to the S-matrix, has the disadvantage 
of not being gauge invariant. A possible solution to this problem is based on 
the introduction of a background field DEI El El. The resulting modified ef¬ 
fective action is then gauge invariant with respect to gauge transformations 
of the background field. This so-called background held method was invented 
to simplify theoretical considerations and also quantum computations in gauge 
and gravitational theories. Originally, the formalism was only applicable up 
to one-loop order. The extension to higher loops was found in DEIEIIHI A 
proof of the renormalizability of Yang-Mills theories in this formalism was given 
in Eiiini, whereas [niiniiisi demonstrate that one obtains the same S-matrix 
as in the usual formalism. The well-known renormalization group functions of 
an arbitrary nonabelian gauge theory at the one-loop [OEIIIIII and two-loop 
level OHIElinilEIilEI] were reobtained in a simpler manner in the background 
held formalism For calculations in quantum grav¬ 

itational theories and nonlinear sigma models, the background held method 
is indispensable. By means of this method Einstein gravity was shown to be 
one-loop m but not two-loop ESOOIEIl hnite. Recently, the background held 
method has found application in the standard model as well ESESIEIISSIESI 

The main purpose of this paper is to demonstrate the effectiveness of the co¬ 
variant background held method in the explicit calculation of the three-loop 
/3-function for a pure Yang-Mills theory. The successful completion of this cal¬ 
culation gives us hope that the open issue of the three-loop renormalizability 
of supergravity may be answered with similar methods. The three-loop renor¬ 
malization group functions for QCD were hrst obtained with conventional held 
theoretical methods in EZIEHI Very recently these results were recovered and 
extended to include scalar helds and Yukawa couplings in m- Although these 
authors use the background held method, they rely upon momentum space 
methods which violate background gauge invariance at the intermediate level. 
Indeed, although the background held method leads to a gauge invariant effec¬ 
tive action, it is common practice to give up background held gauge invariance 
at intermediate levels of a computation □El The reason for this can be traced 
to the desire to treat the interaction with the background held as a perturba¬ 
tion. One commonly chooses the background held covariant Feynman gauge, 
but in order to be able to dehne the gluon propagator and the Feynman rules 
one is led to split up each covariant derivative into an ordinary derivative plus 
a background gauge held term. As a consequence the background held appears 
explicitly in the Feynman rules and manifest background held gauge invariance 
is lost in the process. In addition this splitting up of covariant derivatives leads 
to a considerable increase in the number of vertices as compared with their 
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number in the absence of a background field. Nevertheless, a major gain is 
found in the fact that it now suffices to calculate the two-point function of the 
background gauge field in order to determine the charge renormalization and 
thus obtain the /3-function with relative ease. The contributing Feynman graphs 
are at most quadratically divergent on dimensional grounds (see fig la). At the 
end of the calculation one adds up the various contributions and should obtain 
a transverse answer, i.e. proportional to the square of the background gauge 
field strength. This is usually seen as a valuable test on the correctness of the 
calculations. However, the above procedure contradicts the spirit of the back¬ 
ground field method and leads to unnecessary labor. This becomes especially 
important at higher loops in quantum gravitational theories. 



Fig. a Fig. b 


Figure 1: Comparison between the ordinary (a) and the covariant (b) back¬ 
ground field method. 

In this paper we shall therefore insist on manifest background field invariance 
and avoid the splitting up of background covariant derivatives. This prohibits 
the otherwise standard transition to momentum space and as a consequence all 
our calculations will be performed in configuration space. Such a procedure has 
been advocated before in finding the finite part of the effective action for a pure 
Yang-Mills theory in a covariantly constant background field at one-loop HOI 
ED 1321 and at two-loop m order. It was also used in E2 to determine the 
two-loop /3-function for SQED in a very efficient way. The background field 
will appear only through covariant derivatives in the gauge-fixed action and 
Feynman rules. The number of vertices is thus the same as in the absence of 
the background field. We shall work in the background covariant Feynman gauge 
and obtain background field dependent gluon and ghost propagators via a heat 
kernel representation. Out of these exact propagators and the vertices we may 
construct a compact though formal expression for the effective action at any 
particular loop order. In the chosen gauge there exists a Ward identity |45l I4ti| 
connecting the exact gluon propagator G^v{x,x') and exact ghost propagator 
G{x,x'), namely 

D^G^i,{x, x') + Di,>G{x, x') = 0 

We will use this identity to convert ghost graphs into gluon graphs in the for¬ 
mal effective action and thus roughly half the number of graphs. In order for 
the exact propagators to transform gauge covariantly at both endpoints, they 
contain in particular a Schwinger phasefactor 


^{x, x') 


P exp 
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where the symbol P denotes path ordering. As this phasefactor is dimensionless 
it would seem to lead to quartically divergent graphs. This can be avoided by 
choosing the so-called Fock-Schwinger gauge x^^B^{x) = 0 (note that one is 
free to select different gauges for the background and quantum fields). It is 
well-known that in this gauge the gauge held can be expressed in terms of 
its own held strength. It was shown in rniT^ that one may obtain in the 
Fock-Schwinger gauge an explicit expansion of the propagator in powers of the 
background held strength and its covariant derivatives. In fact, if one demands 
the background held strength to be covariantly constant, then there exists a 
closed expression for the associated heat kernel. The condition DpFfj,^, = 0 is 
analogous to Schwinger’s choice of a constant electromagnetic held strength for 
the electron EiEni With these choices the background held appears in the 
Feynman rules only through its held strength and hence we obtain at most 
logarithmically divergent Feynman graphs (see hg lb). Here, each individual 
graph is background gauge invariant. This is a major simplihcation compared 
to the method described above where one hnds quadratically divergent graphs 
and only the sum of all graphs is gauge invariant. 

To regularize the Feynman graphs we use dimensional regularization with (mod- 
ihed) minimal subtraction [SHESEnilSI Our renormalization procedure is 
based on the i?*-method [SHI Eni' This is a generalization of the well-known 
i?-method EZI which not only eliminates all UV subdivergences, but also all 
IR-divergences for any Feynman graph. Thus the i?*-method allows one to 
renormalize each graph separately without introducing explicit counterterms. 
It can be shown that the /3-function is the same in all MS schemes based on di¬ 
mensional regularization to any loop order m We note that this also holds at 
the level of individual renormalized graphs and this fact constitutes a valuable 
check on our calculations. 

All calculations in this paper were performed with Mathematica isni. At first 
sight, this would not seem to be the obvious choice as Mathematica is not op¬ 
timized for speed or very large expressions, its strength being in the domain of 
pattern matching. However, with the necessary precautions, this programming 
language can handle several thousand terms and finish a calculation in a reason¬ 
able amount of time. Our manifestly covariant approach prevents one from ever 
having to deal with more than a two-thousand terms and hence Mathematica 
works well. 

The paper is organized as follows. In section (21 we give a brief introduction of 
the background field method applied to the Yang-Mills theory. Assuming that 
the field strength is covariantly constant and that the quantum held is given 
in the Feynman gauge respectively the background held in the Fock-Schwinger 
gauge, the corresponding heat kernels and propagators are deduced in sectionjHl 
Pursuant to our covariant approach the resulting Feynman rules are specihed 
in section ^ In section (HI the one-loop calculation is presented. It shows up 
that for determining the UV-divergences it suffices to know the propagators. 
In section IHl we demonstrate that the whole two-loop calculation can be done 
by evaluating only one Feynman graph. Eventually we outline in sectional our 
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three-loop calculation done with the help of Mathematica. Finally we give in 
section!^ an outlook how e. g. a three-loop calculation for supergravity may be 
accomplished by using the proposed covariant background field method. 


2 The background field method in Yang-Mills 
theory 

We consider a nonabelian gauge theory with gauge fields 
ta being the generators of a semisimple gauge group G. The covariant deriva¬ 
tives + Afj, transform homogeneously under a local gauge transfor¬ 
mation U{x), i.e. U~^Dfj_U. The same holds for the field strength 

= [Dfj^ , D^] which appears in the action 

where g is the gauge coupling constant and we normalized the generators in the 
representation R by tr{tatb) = —T{R)Sab- We shall work in euclidean space 
throughout. The equations of motion and Bianchi identities are given by 

= 0 , DpF^i, + DpFyp F D^Fp^ = 0 ( 2 ) 

We split the field into a background field Bp and a quantum field Qp such 
that 

Ap = Bp F Qp ( 3 ) 

This induces a similar splitting for the field strength 

= Fpi, F DpQi, — DjjQp F Qv\ (4) 

where on the right hand side Dp and Fpi, represent the background covariant 
derivative and field strength, respectively. From here onward this notation will 
be left understood. An infinitesimal gauge transformation of Ap with parameter 
A can of course be distributed in many ways over Bp and Qp, but the two most 
convenient choices are the “quantum transformation” 

6Bp=0 , SQp = DpAF[Qp,A] (5) 

and “background transformation” 

dBp= Dp A, SQp= [Q^,A] . (6) 

The trick is to add now a gauge-fixing term which breaks the quantum gauge 
invariance, but respects the background gauge invariance, i. e. we add to the 
classical action a term 

*5fix = ^ / dx (7) 
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where the gauge-fixing function G transforms covariantly under The asso¬ 
ciated Faddeev-Popov term is 


/ SG 

dxtvC*—C , (8) 

where the variation of G under a quantum gauge transformation is meant and C 
and G* are the ghost and antighost fields. The total action ^tot = <S'ci-l-S'fix+<5'FP 
then appears in the generating functional for all Green functions 




N 


VQ VG* VG exp(-S'tot + J ■ Q + v* ■ G + G* ■ r^) (9) 


with sources J, t] and t]* and normalization factor N chosen such 
that ^[5,0,0,0] = 1. Note that the background field is not coupled to the 
source J. The generating functional for the connected graphs is then given by 


W[B,J,r,*,r^]=\uZ[B,J,7^*,r^] . 

Defining expectation values for all fields as follows 

^_SW 

^ SJ ’ Srj* ’ Sr] 


( 10 ) 


( 11 ) 


one defines the effective action by 

r[B,Q,C*,C]= J-Q + r^-C* +C-r]* - W[B,J,r]*,r]] (12) 


The functionals Z and W will be invariant under background gauge transfor¬ 
mations, if all sources and ghost fields transform in the same way as Q does. 
The same holds for T if one demands Q, C and C* to transform as Q. A good 
gauge choice is e. g. the generalized Gervais-Neveu gauge 

G = , (13) 

but we will restrict to the background covariant Feynman gauge a = l,/3 = 0 
with associated Faddeev-Popov term 


5fp = 2 j dxtv{G*D^^{D^G+[Q^,C])) . 

If we now, as mentioned briefly in m decompose the current in 

Jfi = 3)1+ 3)1 


with 




5Stot[B,Q] 

SQf^ 


Q = 0 


(14) 


(15) 

(16) 
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we can write the total action as 


Stot[B,Q]-J-Q 

= Sci[B] + ^Jdxtr iQf^A^,[B]Q’')- J dxtT{C*A[B]C) 

-S,nt[B,Q]-J-Q (17) 

with the definitions 

^^B] = -{S^,D^ + 2F^,) (18) 

A[B] = -D^ (19) 

S^nt[B,Q] = J dxtY(^{D^Q,)[Qf^,Q''] + ^[Q^,Q,][Qi^,Q''] 

+ C*D^[Q>^,C]) ( 20 ) 

where the field strength is given in the adjoint representation. Following the 
definition of apDendixOthis means that corresponds to the commutator 

[F^u.Q^]- 

Thus we get a new generating functional 


Z[B,J,r]*,r]\ 


j VQVC*VC exp 


-Sci[B]-\ J dxQf^A^^Q'^ 


dxC*AC 


~ Sint [B,Q] + J ■ Q T]* ■ C C* ■ rj 


( 21 ) 


depending on a new current J. The exponent is now free of any term linear in 
the gauge field Q, except the one coupling to the current J. 

Imposing the constraints Q = C = C* = 0, the effective action of m now 
reads 

r[i3,0,0,0] = -W[B,J,r] (22) 

and is still invariant regarding the gauge transformations of the background field 
given in ©■ The pertubative expansion of this effective action ijSl contains 
only vacuum graphs, which are of course IPI. Thus, the number of emerging 
Feynman graphs is drastically reduced as we will show later on. 

The expansion of the effective action in terms of h is given by 

r[B, 0,0,0] = Sci[B]+ (^^lndetAf,4B]-lndetA[B]^h + 0{fi‘^) . (23) 

Hence the effective action r[i3,0,0,0] equals the classical action Sci[B], if we 
neglect all quantum effects by using the approximation ?i —> 0. To calculate the 
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correction linear in Ti which corresponds to the ordinary one-loop calculation, 
we need to evaluate the logarithm of the determinant of the operators A^i, and 
A. Higher corrections in fi are represented by Feynman graphs, which are built 
according to the Feynman rules of sectional 


As expected while expanding the new effective action r[i3,0,0,0] we meet the 
problem of UV-divergences. Following the prescriptions of the background field 
theory it suffices to redefine the background field B and the coupling constant 
9 by 


^ B^^ 

9 ^ 9o = 


(24) 


The quantum field Q and the ghost fields c* and c however need not to be 
renormalized, since they can only be found inside of the vacuum graphs and as 
a result the renormalization constants of the vertices would be cancelled with 
those of the propagators |7|. 

The renormalized field strength is given by 

= Z]^^(d^B,-d,B^ + zF^Z]^^[B^,B,]j, (25) 

and has to be invariant in terms of the background field transformation 
Thus the two renormalization constants have to satisfy the identity 

Zi/24/2 = 1. ( 26 ) 


As already mentioned in this is the reason why we have only to evaluate the 
two-point functions of the background held in order to renormalize the Yang- 
Mills theory. Preserving the covariance, as done in the presented calculations, 
leads directly to two-point functions of the held strength and has moreover the 
advantage that, in comparison to the common procedure, the number of vertices 
does not increase. The corresponding Feynman rules, given in section^] contain 
two 3-vertices and one 4-vertex, whereas the Feynman rules of m include four 
3-vertices and hve 4-vertices. Hence the preservation of the covariance of the 
background held reduces the number of vertices and consequently the number 
of Feynman graphs. 

In the case of the two-loop calculation the number of Feynman graphs can be 
reduced once more from twelve in |22| to one as shown in sectional 

A nice proof of the renormalizability of the Yang-Mills theory in the background 
held formalism can be found in m Based on the fact that for a vanishing back¬ 
ground held we hnally get the original Yang-Mills theory whose renormalizabil¬ 
ity has already been proven, the authors Liischer and Weisz show, by using the 
BRS-, the background held- and a shift-symmetry, that the Yang-Mills theory 
written in the background held formalism is renormalizable as well. 
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3 The exact propagators 


For a given operator A the corresponding Green function or propagator in d di¬ 
mensions is defined by the wave equation 

AG{x,y) = S{x -y). (27) 

The associated heat kernel K satisfies the heat kernel equation 

+^ K{x,y]T) = {) (28) 

with the boundary condition 

K{x,y,{)) = 5{x -y). (29) 


The Green function G is connected to its heat kernel K via the equation 

poo 

G{x,y)= / dTK{x,y;T) (30) 

which therefore satisfies (|77|i. 

For the ordinary d’Alembertian Aq = —d^ it is straightforward to verify that 

Koix,y;T) = ^_^exp(^-^(a;-2/)2^ , (31) 

which depends only on the coordinate difference x — y, solves (1^ . Thus, gen¬ 
eralizing this result, one makes for the heat kernel associated to the scalar 
operator A of equation m the following ansatz 

K{x,y\T) = 1 /)L(t) exp - y)^M^^(r)(a: - yY^ , (32) 

in which the phasefactor ^{x,y) is defined by 


^{x,y) = Pexp 
Imposing the constraint 


B^{z)dZfi, 


(33) 


implies 


(DaF^^) = [Da,Ff^„] = 0 


[Fpa, Fp^] = 0 . 

Furthermore, Shore shows in m. that the constraint 
identities 


(34) 

(35) 

implies the following 


= Ffj,^{x)<^{x,y) 

= -^Fp^{x){x-y)''^ix,y) 

= -^^ix,y)Fp^{y){x-y)''. 


‘^{x,y)Fp^{y) 

Dpix)^{x,y) 

^{x,y)Dp{y) 
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(36) 

(37) 

(38) 




Inserting the ansatz into the heat kernel equation (UHl and taking into 
account the condition in and the constraint m, one obtains as solution 

exp vTF^aT (cothFr)'"^ {x - yY^ (39) 

in accordance with EH, where Tr denotes the trace over the Lorentz indices. 
The fraction following the Tr is a formal expression, which is defined by its series 
expansion 




\ sinh Ft 


oo 9 (r) 2 k-l _ 1 \ 




{2k)\ 


(40) 


where F^)) is a short notation for the expression F^F^')'^ ■ • ■ A detailed 

explanation of our notation is given in appendix O 

Due to the identities (E3 and EEli the heat kernel of the vector operator 
can be deduced from the heat kernel of the scalar operator A. To satisfy the 
heat equation 

- 2F^a) K\ = 0 

and the corresponding boundary condition 


Kij,^{x,y,0) = I6^^d{x-y) 

it is sufficient to extend the heat kernel of A by an exponential factor. In detail 
we get 


K^u{x,y,T) = (exp2F(a;)r)^^ A:(x,y;r) 

= A:(a;,?/;r) (exp2F(2/)T)^^ , (41) 

where, depending on whether we multiply the additional exponential factor from 
the left or from the right, the field strength in the exponential factor depends 
on X or y. 

The two heat kernels are connected by the following identity 

D>^K^, + KD, = 0, (42) 

sometimes called Ward identity M This identity holds even without using the 
constraint as shown in |4til Iti2| . 

Now, by using PHl . we can determine the propagators corresponding to these 
heat kernels. They satisfy the following heat equations 


AG{x,y) 

^p.\GY{x,y) 
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lS{x-y) 
lSf_,„6{x - y) 


(43) 

(44) 






and are connected by the identity 

D^^G^, + GD, = 0 (45) 

which can be deduced from 63 by integrating over r. This identity will be 
of great help, since all ghost graphs contain scalar propagators with a covari¬ 
ant derivative acting on it. By transforming these into vector propagators the 
number of graphs can be reduced about 50%. 

Applying further a covariant derivative from the right to the identity lUSl) , we 
obtain 

D^Gfj,^{x,y)D'' = -G{x,y)D‘^ = -D^G{x,y) = lS{x-y), (46) 


which is another helpful identity, since it leads to the cancellation of an exact 
ghost propagator. 

Unfortunately, it is not possible to find a closed representation for the propaga¬ 
tors. But as long as we are only interested in the renormalization of the pure 
Yang-Mills theory it is quite enough to take the expansion up to quadratic order 
in the field strength. For convenience we set z = x — y and get the following 
series expansion for the scalar propagator 


G{x,y) = 


dTK{x,y;T) 


^{x,y) J dr Ko{z;t) - 


— Z>^YlX - — F2 
12 12 


0(F3) 


= <^{x,y)G^{z)-^z>^¥l,z''G^{z)-^¥^G2{z) + 0{¥% (47) 

where Go represents the free propagator in coordinate space. By convoluting i 
functions Gq one after another we get the function Gi multiplied by a factor i\ 
(see appendix 0. The vector propagator is given by 

poo 

G^^ix,y) = / dTK{x,y;T){exp2FT) 

Jo 

poo 

= / dTK{x,y,T){S^, + 2F^,T + 2Fiy)+0{F^) 

Jo 

poo 

= / dT[K(x,y;T)Sf,^ 

Jo 

+ Ko{z; r) (2 ^x, y)F^,r + 2 + 0(F3) 

= 5^,G(a;,y) + 2$(x,2/)F^,Gi(z) + 2F2 G2 (z) + 0(F3) .(48) 


The phasefactor $ depends on the background field B and for this reason de¬ 
pends on the field strength F. Hence the phasefactor can be set to one in all 
terms containing the field strength in quadratic order. 
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Using the identities in appendix m the explicit coordinate difference can be 
replaced by derivatives of the Green functions. Doing so, the series expansion 
of the scalar propagator is given by 

G{x,y) = ^{x,y)Go{x-y)-^F^G 2 ix-y) 

-^Fl,d^d^G3{x-y) + OiF^), (49) 

respectively in graphical notation by 

- = ^x,y), -. . . 

-^F^,6. > I - + 0(F3). (50) 

The series expansion of the vector propagator is given by 

G^j.v{x,y) = i^{x,y)5f,^Go{x - y)+ 2i^{x,y)F^^Gi{x - y) 

^}-^{SF%-FH^.)G2{x-y) 

-iFlJ^,d^d^G^{x-y) + 0{F^), (51) 

whereas the graphical notation is of the shape 

- - - = <^{x,y)5f,^. - • + 2^{x,y)Ff,^. -.-. 

+ 1(8F^,-FV)2. ■ ■ ■ 

- g- + 0 (f3) . (52) 

As shown above the graphical representation of the exact scalar propagator is 
illustrated by a thick gray line and the exact vector propagator by a thick black 
line. Unfortunately the remaining phasefactors prevent our expansions of the 
propagators to be translation invariant. If, however, we use a procedure given 
in m we can restore the translation invariance for the whole Feynman graph 
as shown in chapter El and 0 

While analyzing the expansions (133 and m of the propagators, we realize 
that any possible divergences must reside in the Green functions. Following the 
argumentation of appendix 0] the functions Gi contain no UV-divergences but 
for d = 4 — e the functions Gi with i >1 have to be replaced by the functions 
Ri in order to cancel the IR-divergences. As result in 4 dimensions the function 
i?i is the only one which harbors an IR pole with constant residue. Therefore 
on the diagonal the scalar propagator is given by 

G{x,x) = 0 (53) 
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whereas the vector propagator satishes on the diagonal the equation 


which will be needed to deduce the value of the tadpole defined in section 0 

After imposing the Feynman gauge (a = 1) and the constraint (IHdll we are still 
free to choose a gauge for the background held B^{x). Instead of using the 
Feynman gauge we select the Fock-Schwinger gauge I^Sl 

x>^Bl{x) = 0. (55) 

Later on we will see that this gauge allows us to transform every generated 
Feynman graph into logarithmic divergent Feynman graphs. As a result the 
divergences of these graphs may be represented by a Laurent series in coordinate 
space containing no derivatives, respectively in momentum space containing no 
dependence on any momentum. 

As an aside, we note that by applying a Fourier transformation to the Fock- 
Schwinger gauge we are lead to 

d^Bl{p) = Q . 


In this sense the Fock-Schwinger gauge in coordinate space and the Lorentz 
gauge in momentum space are dual to each other. 

As a result of using the Fock-Schwinger gauge, the background held can be 
represented as a function of its own held strength [23I1HI To prove this, we 
take the derivative of the Fock-Schwinger gauge 


B„{x) + x^d„Bf_,{x) = 0 

and replace afterwards 9,^5^ by —-I- d^B^, + 5^]. Since by dehnition of 

the Fock-Schwinger gauge the expression x^[B^{x), Bij^{x)] vanishes, we hnally 
get 

(1 -h x^d^)B^{x) = x>^F^^{x). 

Using the transformation x = Xx and the identity Xd/dXf{Xx) = x^9^/(Ax) 
we are led to 

+ S,(Ax) = Xxf^F^^iXx) 

with solution ^ 

B^{x) = [ dWx^Ffj^^{\x) (56) 

Jo 

Since on the one hand the background held in the Fock-Schwinger gauge is 
dehned by its own held strength and on the other hand the commutator of two 
held strength tensors vanishes, the constraint m reduces to 


dpF^„ = 0 , 
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which means that the field strength is constant. Finally this allows us to perform 
the integration in equation and we get a potential which depends linearly 
on the field strength 

B,{x) = (57) 

from which we deduce the following phasefactor 

<^{x,y) = exp (^x^F^^y''^ , (58) 

which does of course satisfy the equations 18tii:-{8ll . Later on we will see that 
this special form of the phasefactor is an important requirement for the trans¬ 
formation of all occuring Feynman graphs into graphs which are at most loga¬ 
rithmically divergent. 


4 The vertices 

After the transformation Q —> g Q the Feynman rules for the Yang-Mills theory 
are given by the following expressions. 

The triple gluon vertex can be written as 

g J dx ( f‘^"^ 6 prS{x-X 2 )Six-X 3 )D^‘^S{x-Xi) 

- f^‘^'^5pa5{x — X2)5 {x — X3)D'^°‘5{x — Xl) 

+ r'^5^r5{x - X3)5 {x - X 3 )D^H{x - X 2 ) 

- r^5p^5{x - X3)5 {x - X3)D'^^6{x - X 2 ) 

+ r^-^S^rSix - xi)d{x - X 2 )D'^‘^ 5 {x - X 3 ) 

- r^'^5pr5{x - X3)5 {x - X2)DZ5{x - X 3 ) ) 


p, a cr, 6 



T,C 
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The arrow indicates a covariant derivative operating on the line on which it is 
drawn and the circular arc denotes the contraction of its Lorentz index. As the 
lorentzian indices of the two remaining propagators are always contracted the 
second circular arc can be omitted. 

Concerning the triple vertex containing two ghost fields, we have to extend our 
graphical representation 


a b 

/x,c 


(60) 


by an additional rule which says, that only those graphs are taken into account 
where each gray line with a derivative is connected to a gray line without. 

The quadruple vertex is dehned by 

f dx S(x — xi)S(x — X 2 )S(x — X 3 )S(x — X 4 ) X 



T,d a,c 


Note, that regarding the graphical representation, we have inserted a i5-function 
represented by a dashed line in the quadruple vertex. The two circular arcs 
indicate, which indices of the two triple vertices are contracted. Again the two 
remaining indices are contracted as well. In those cases, where the two gluon 
lines of each vertex are contracted, we omit the arcs entirely. 
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Thus, as an advantage of this representation, graphs containing only triple ver¬ 
tices are generated, which vastly reduces the amount of possible topologies. 


5 One-loop 

It is well-known that the divergences in the order of one-loop can be determined 
by the heat kernels (see i. e. ESI)- In the case of pure Yang-Mills theory we are 
led to 


hi 


1 det det A[.B] 

2 detA^,y[0] detA[0] 


lim ds 
s^O 



dr' 


tr 


dx 


x; t\B) - K^^^{x, X] t|0)) 
— 2(^K{x, x]t\B) — K{x,x\t\Q)) 

2(^1 * 71 ^ 



(62) 


As a matter of fact this integral cannot be solved. But since we are only inter¬ 
ested in the UV-divergences we expand it in terms of r and after the following 
integrations take the limit r ^ 0. 

Now with the help of 


exp 


1 j / sinhFr 
~2 Fr 




1 

12 ' 


1 -k —FV^ -k 


and 

(exp2Fr)^^ = dl-2F^T^ + ■■■ 

we eventually get 

which corresponds to the results of the literature nmsi. 


(63) 

(64) 

(65) 
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6 Two-loop 


The contribution of the quadruple vertex (ED to the two-loop effective action 
is given by 




+ 




( 66 ) 


.0 _ 0 


By using the identity lE3J), it can be shown, that the second graph in the first 
row is similar to the Hrst graph in the first row. Due to the same identity llMll 
the expression x) yields null and hence the first graph in the second row 

and the second graph in the third row vanish. If we remove the twist in the 
first graph of the second row, we get an additional minus sign as a result of 
the antisymmetry of the structure constant. Finally the two remaining graphs 



Figure 2: Tadpole identity 


are shown to be equal by means of the tadpole-identity illustrated in figure |21 
Written as a formula we get 


fabcGf^(x,x)S(x 


x')fab'c'G^,,:ix\y')G'^,,:{x',z') = 

2G%{x,x')6{x - x')fabcGr,i^,y')fac'b'G'',,:{x',z') (67) 
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which can be proved by using the Jacobi identity and equation (15411 only. 

Thus, at least the contribution of the quadruple vertex shown in (TO) can be 
reduced to a single graph 



where the minus sign on the right hand side is a result of equation IS3J). All 
graphs sharing the topology of those in equation inn) will, according to their 
shape, from now on be called 0 graphs. 

Following (1^ the contribution to the two-loop effective action of the two triple 
gluon vertices is given by 36 0 graphs. Fortunately the symmetry of the 0 graph 
allows us to replace the sum (ICTll of one vertex by one of its terms multiplied 
by 6 without loss of generality, leaving us with only 6 graphs. Furthermore 
one of these graphs can - by exchange of the two coordinate space points - be 
transformed to a graph already existing. Adding hnally the contribution of the 
ghost vertices (IHHll we get the two-loop effective action as shown in hgure El 
Thus, using the tadpole identity and some symmetries we are left solely with 
graphs of the 0 topologies. This simplihcation gets even more important in 
higher loop calculations, since the amount of possible reductions is increasing. 

Before proceeding we will establish a notation, which will be just as powerful as 
the graphical representation. As very useful we consider the bracket notation 
introduced in EH. This notation is free of group indices and the symmetries 
are easily seen. It is defined by 


A 





(69) 


c 


where - due to clockwise representation of the two structure constants - the 
graphical representation holds a minus sign. 


Thus, the symmetries are given by 


(A,B,C) = (A,C,B) = (B,A,C) , (70) 


which means that the three propagators of a 0 graph are symmetric to each 
other. 

The bracket notation corresponding to the graphical representation of the effec- 
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Figure 3: Two-loop effective action 


tive action of figure 0 is given by 
3 


r, = 


dx dy 


(1; G/il/, 


+ 2 p<7•) DpGfxif Da) + {DpGfjv^ Gpa, GpvDcr) 

+ {Gpu, DpGpvD a, Gpa) — 2 {GpvD rj, DpGpu, Gpa) 
+ {GpuDa, Gpu-, DpGpa) 


~ o (DpG, Gpv, GDy) 


(71) 


One should bear in mind, that up to now only the Feynman gauge has been 
used. 
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But the effective action GD can even be more simplified. In a first step we 
transform with help of the socalled Ward Identity (EH) all ghost propagators 
into gluon propagators. 

The graphical and the bracket notation are illustrated by 



^ {D^G,G^,,GD,) = {G^,D,,G^,,DpGp,). (72) 

The new emerged graph has a remarkable feature: All covariant derivatives act 
on the propagator with which their lorentzian index is contracted. Such graphs 
will from now on be called primary. Concerning 0 graphs with one derivative at 
each vertex there exists still another primary graph. It is generated by replacing 
the 5-function according to equation EU in the first graph of figured Thus we 
get 



^ = {Gp,,DpGp,D,,Gp,). (73) 

Before introducing further simplifications we have to show that the covariant 
derivatives follow the rules of partial integration. To prove this statement we 
take the sum of three vertices resulting from a partial integration and get 

{DpGi,G 2 ,G 3 ) + {Gi,DpG 2 ,G 3 ) + (Gi, G2, Z^^Gg) 

= dp(Gi,G 2 , G3) -I- {[Bp, Gi], G2, G3) -I- (Gi, [Bp, G2], G3) -I- (Gi, G2, G3]) 

= 9^(Gi,G2,G3) + [i?^,(Gi,G2,G3)], 

where the last equal sign stems from the identity D[G, E] + [G, D]E = [G, DE], 
Eventually this result yields null, since the commutator vanishes due to the fact, 
that the expression (Gi, G2, G3) is a scalar and the total derivative vanishes due 
to boundary conditions. 

As shown in figure 0] by partially integrating each vertex of the two primary 
graphs we get the remaining five graphs with correct coefficient, but sometimes 
incorrect sign. 

The correct sign can be deduced from the number of arcs which connect the 
derivative and the associated propagator clockwise. For an odd number the 


20 






2 



Figure 4: Partial integration applied to the primary graphs 

sign is negative. This procedure in what follows will be called star operation 
and will be indicated by a symbol at the right upper side of the graph. The 
graphical representation of the star operation as applied to both primary 0 
graphs is illustrated in figure 0 



Figure 5: Star operation applied to the primary graphs 

Thus the 0 graphs of the Yang-Mills theory can be divided in two classes. By 
partially integrating the two vertices of the primary graph with both derivatives 
operating on the same propagator, we generate all graphs belonging to the first 
class. The associated graphical representation can be found in the first row of 
figure 0 The second primary graph gained from the ghost graph by means of 
(d, leads us by partial intergration to the graphs of the second class (as shown 
in the second row of figure E)). 
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Eventually the two-loop effective action by means of the star operation, and the 
identities lB51) and (EH) can be reduced from originally seven graphs down to 
three as illustrated in figure El 

It should be emphasized that up to this moment we have only used the Feynman 
gauge since, as mentioned before, the two propagator identities (1^ and 
are even valid without the assumption DpF^i, = 0. 

It should be mentioned that the gauge invariance of the background field 
is preserved, and since all identities used until now are non-pertubative as far 
as the background field is considered, the new effective action can be used to 
calculate the two-loop finite part as well. 



Figure 6: Simplified two-loop effective action 


Assuming in addition that the field strength is covariantly constant m we are 
now allowed to use the expansion of the propagators EH EH). All terms of 
higher than quadratic order in the field strength are dropped. 
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After considering the effect of the covariant derivative on the phasefactor (|!r7 | |!TH |i 
the latter can be substituted by the identity in all terms containing the field 
strength quadratically. 

Imposing the Fock-Schwinger gauge on the background field we get the 
phasefactor defined in and a field strength which is constant. 

Due to the translation invariance we are allowed to fix a vertex without restric¬ 
tion in general at the origin of coordinate space gHlEll. Thus all remaining 
phasefactors in the two-loop graphs can be neglected, meaning 

$(0,2/) = 1. (74) 

Now all phasefactors have disappeared and we see immediately that graphs 
containing the field strength linearly do, for symmetry reasons, not contribute 
to the result. Graphs without field strength are neglected. 

Finally only graphs provided with the following properties are left: They contain 
no phasefactor at all, but the field strength in quadratic order, and they are log¬ 
arithmically divergent. Since the phasefactor has been the only term depending 
on points in coordinate space as well as on some group indices, the remaining 
graphs can be decomposed in an ordinary Feynman graph and a graph consist¬ 
ing of structure constants. A step by step calculation of the second and third 
graph in the last row of figure IHl is given in appendix ini However a sketch of 
this procedure applied to the last two graphs of figure El is given by 

{G^v,GpvDrj,DpG^a) = (‘hGo, $GoD^, iD^$Go)-l-■ • ■ 

= -($Go,$a^Go,$a^Go) 

+ (‘I’Go, $F^pcIpGi, ^Fpo-^o-Gi) -!-••■ 

(Go,FppapGi,Fp,,a,Gi) + --- (75) 

The last equal sign stems from ( 171 ) . From the 16 graphs which we get after the 
expansion of the propagators, we will calculate as an example the graph which 
depends on the field strength only via the covariant derivative and the phase 
factor The remaining graphs are transformed to logarithmically divergent 
graphs containing the field strength quadratically in a similar way. 

(Go,FppapGi,Fp,,a,Gi) = Goix-y)dpGiix-y)d,Giix-y){l,Fpp,Fp^) 



(76) 
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The first graph represents just an ordinary Feynman graph. The corresponding 
divergences can be looked up in the table of appendix^ whereas the value of 
the second graph consisting of structure constants can be deduced by means of 
the rules as listed in table |0 After determining the ordinary graph, the lacking 
contractions have to be made. The following two identities 


(1,1, F^) = 

G2tr(F2) = 

gI{.f;S 

(1,F^,„F^^) = 

iG2tr(F2) = 



suffice to obtain the result for all remaining two-loop graphs made of structure 
constants. 

Resisting the temptation to calculate the Feynman graphs immediately, we show 
that by using several manipulations given in appendix^ the number of graphs 
can be reduced to two and finally even to one. 

Thus the two-loop divergent part of the effective action is given by 

pdiv 
^ 2 


1 ) Gfiv) — {Gfj^v, Gpa, D pG pa D n) — 2{G pin Gp,yDa, D pG pa) 

- ^(Go,Go,GiF 2) + 0(e°) (77) 

11 D 



= ^ . (1,1,F2)- ^ -. (l,l,F2) + 0(e°)(78) 

= -^(Go,5;.Gi,5'^GiF") + 0(e°) (79) 


^ (1,1,f 2)-pO(e°) (80) 


whereas the transformation from the graphical representation iTTHll to re¬ 
spectively from equation (EZI) to GHJ is accomplished by using the identity itTT^ . 
The resulting graph CZHSOI) has a remarkable property: it contains no subdi¬ 
vergences and therefore has only a simple pole. Thus just by looking at the 
hnal graph without calculating any integral, we realize, that the quadratic pole 
of the Yang-Mills theory vanishes at the order of two-loop. 


24 





7 Three-loop 


In the following we will generalize the methods developed in the last section 
such that they are suited to determine divergences of Feynman graphs beyond 
two-loop. Afterwards we will by means of these methods calculate the /3-function 
of the pure Yang-Mills theory up to three-loop. 

After using the tadpole identity (see figure|2J in three-loop the remaining graphs 
belong to one of the following two topologies. The graphs of the first topology 
will be called ball graphs, because of their shape. The bracket notation, respec¬ 
tively the graphical representation is given by: 


{A,B,C\K,L,M)x 


fa ibiCi fa2b2 C2/a-3^3C3 fa4b 


4C4 


j^aia2 ^b 2 63 (JC2C3 04 J^b 3 64 j^C3C4 



(82) 


The ball graphs obey the following symmetries 


{A,B,C\K,L,M)i 


flip vertical 
flip horizontal 


{A,C,B\K,L,M)i 

(k,b,c\A,L,m), 

iA,L,M\K,B,C)u 


whereas propagators with a tilde have to be read backwards {G^v{x,y) = 
x)). Thus we are left with 16 possible representations for each ball graph. 

Due to their symmetries the graphs of the second topology are called tetrahe¬ 
dron graphs. We introduce the following bracket notation respectively graphical 
representation: 


{A,B,G\K,L,M)2 


faibici fa2b2C2 fa^b^c^ f 

j^aia 2 j^bi &4 ^ciC3 j^atas ^6362 j^C 2 Ci 
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Since the tetrahedron graphs show up the following symmetries 


{A, B, C\K, L, M)2 {B, C, A\L, M, K)2 (C, B, A\M, L, K)2 

rotated M, i|B, C, i^) 2 , 


each graph has 24 representations. 

After using the Ward identity to eliminate all ghost propagators the three- 
loop effective action consists of the ball and tetrahedron graphs found in figure[71 
and|5| 
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(' 3 ') 

The fifth graph in the third row of is a hybrid between a ghost graph and a 
star graph. As shown in figure El in the original graph the lower loop is built up 
only by two ghost propagators, whereas the upper loop is built up by two vector 
propagators with the star operator acting on it. A graphical representation and 
the corresponding bracket notation can be found in figure El 






Figure 9: Hybride ball graph 

The second graph in the last row of in hgureElneeds some explana¬ 

tion also. It is built up by a ghost loop containing three ghost propagators and 
a vertex with a star operator acting on it. A graphical representation and the 
corresponding bracket notation can be looked up in figure [Till 
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Figure 10: Hybride tetrahedron graph 


Taking only graphs which contain no (^-function, we see that the ball graphs can 
be divided in eleven classes, whereas the tetrahedron graphs can be divided in 
six classes. 

The graphs can be categorized in accordance with corresponding primary graphs 
as shown in the first two columns of the tables ^ and |3 The third column 
holds the number of graphs we get when applying the star operator at each 
primary graph. By transforming the ghost propagators of the ghost graphs into a 
gluon propagator (E31) we obtain primary graphs which after partial integration, 
significantly reduce the number of graphs in the associated class. Regarding 
classes with no matching ghost graph, we add a suitable primary graph, perform 
the partial integration and eventually substract the same primary graph. Thus 
the number of graphs containing no 4-vertex can by means of this procedure be 
reduced to nearly half the number as quoted in the fourth column. The numbers 
written in brackets signify the amount of artificially added primary graphs. 
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Table 1: Classes of ball graphs 

Thus, as shown in tabled just by using the Ward identity respectively adding 
some well chosen primary graphs, we reduce the number of ball graphs (com¬ 
posed of 3-vertices) from originally 99 to 49. The remaining ball graphs contain 
at least one 4-vertex and therewith at least one (5-function. 

As shown in table [3 by using the Ward identity respectively by adding primary 
graphs, the number of tetrahedron graphs can be also reduced significantly, 
namely from 60 to 31. 

The remaining ball and tetrahedron graphs contain at least one ^-function. 
Therefore, by means of the Jacobi identity such ball graphs can be transformed 
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Table 2: Classes of tetrahedron graphs 
in two tetrahedron graphs as shown graphically in ilMli . 



The value of the graphs consisting of structure constants can easily be deter¬ 
mined by use of table 10 In particular we get: 


(1,1,1|1,1,F2)i 

- CliFp^. 

(1,1,1|F2,1,1)i 

- CUFp^. 

(1,1,1|1,F^^,F^^)i = 


(l,l,l|F^^,l,F^^)i = 

1 /^3 ( ipa \2 
~2 ^2 ■ 

(l,l,F^^ll,l,F^^)i = 

-ici{F-^r. 

(F^^,l,l|F^^,l,l)i = 

- CliFp^: 

(1,1,1|1,1,F2)2 

1 /^3 f TPa \2 

2 ^2 : 

(1,1,1|1,F^^,F^^)2 = 


(1,1,F^^|1,1,F^^)2 = 

0 


(85) 


( 86 ) 
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Following our prescription of the two-loop calculation we expand according to 
the propagators in each graph. After discarding all graphs which depend 
on the field strength more than quadratically, we set the phasefactors equal 
to one in those graphs containing already two field strength tensors. In the 
remaining graphs we let the covariant derivatives take effect on the phasefactors. 
Since the phasefactors depend on points in coordinate space, whereas the 
propagators depend only on differences of these points (translation invariance), 
we will rewrite the phasefactors in the following manner. Take an arbitrary path 
xo, • ■ • ,Xn of a graph connecting the point Xn with xq- Thus the phasefactor 
connecting the last two points of this path may be written as 

^(Xn-l,Xn) = exp 

‘ fT, 

^ \i^0 


^ n—1 

“ ^ ((Xn-2 

^n—l) ) 

^ n—1 

g ^ ^ ((^n—2 ^n—i—1 (^n l) “1“ 

2^0 





3^71—2—1 (^n ^n—l) “1“ 


XQFfj_„{Xn - Xu-iTY + 


(87) 

Following the instructions of |48| respectively our two-loop calculation, we fix xq 
at the origin of our coordinate space by setting its value to zero. As a result the 
phasefactor of equation ilSTI) depends solely on differences of points in coordinate 
space. Using the identities of appendix lEI these differences can be absorbed in 
the Green functions of the propagators. Finally, discarding again all graphs 
not containing the field strength quadratically, we have decomposed the exact 
graphs into ordinary Feynman graphs each multiplied by a product of structure 
constants. 

We want to reemphasize, that both, the fixed point and the path, can be chosen 
freely since the divergences of an exact graph do not depend on the choice of 
these. 

After expanding the exact graphs by means of we get, after using some 

obvious identities, about 2000 graphs. Again, we want to refrain from calculat- 
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ing all these graphs explicitly, but reduce their number by using similar identities 
as demonstrated in the two-loop case. 

First we eliminate the S tensors by help of the identities of appendix^ By doing 
so, we have to watch carefully, if the S tensors, building up a dimensional d, 
belong to a divergent subloop of the graph. 

Thereafter we gather all tetrahedron graphs containing four derivatives. Then 
we build up a basis which contains only graphs without subdivergences, like the 
one given in figure ITTl Since these graphs have only simple poles, we can use 
the identity 

( 88 ) 

while neglecting the e part of the dimension d, which means d equals four. 



Figure 11: Calculation of a simple tetrahedron graph 
In figure im we give the most elementary example of a tetrahedron graph with 
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four derivatives being reduced to the ordinary tetrahedron graph without deriva¬ 
tives. 

Applying a S tensor to equation (ISSt we get the identity 

= (89) 

which allows us to reduce now tetrahedron graphs containing two derivatives, 
but no subdivergences to the ordinary tetrahedron graph without derivatives. 

Thus, after some algebra the sum of all tetrahedron graphs should be simplified 
in such a way, that only the ordinary tetrahedron graph and some other non- 
tetrahedral graphs remain. After doing this calculation we see, that even the 
coefficient of the ordinary tetrahedron graph vanishes, which explains why the 
final answer contains no <('(3) term. 

At present we have to calculate finally 24 Feynman graphs to determine the pole 
part of the effective action. Unfortunately, due to the amount and complexity of 
the rules needed for doing this task automatically, some parts of this reduction 
were introduced by hand. In the following, we give our complete answer for the 
effective action. Each graph is illustrated in a pair of brackets. The first row 
contains the graphical notation where the following definitions are used: 

d/, 

= (1,1,1|1,1,F2J2 
= (1,1,1|1,1,F2)2 . 

The bracket notation is to be found in the second row. The divergences of the 
graph are given in the third row where we omitted a factor (7|/(167r^)^S'ci. The 
definition of this bracket can be found in appendix IXI 



pdiv _ 

^3 — 


/ 


1 

6 


V 


\ 



(Gi,Go,dp^Gol 

1, 9pctG3F^j,)2 

(0, 0 , 1 ^ 2 ) / 


p2 

/i.1^ 

(Go,Go,9pGo| 

1, Go, i9ppyG3F^,^)2 

V (0>0,rk) 
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p2 

fiiy 


( 1 , Go, dpi^Gol 

l,G2,9p.G2F2j2 


V ( 0 ,-ii) J 

( \ 



(l,Go,Go| 

1, Go, 9pi,G3F^^)2 


V ) 


( \ 



(1, Go, Go 

1, dpiyGo, G3F^j,)2 


V (0^0^ m) J 


( \ 



(1, Go, l9p;yGo 
l,Gi,G2F2j2 


/3507715 16281731 , 18422737 ,, 

-fiH- 

\ 5184 414272 27648 



V 
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/1093061 36075, 41537,0 


V 4608 


1024 


9216 


+ - 


1180189 

13824 


242191 , 901981 ,o 

-d - ———fi 


4608 


165888 



(l,Go,Go| 
1,Gi,GiF2)2 

/ 


10009 

27648^ 


(l,Gi,Gi| 

1,1,GiF2)2 

V (- 4 , 0 , 0 ) y 


+0(e 


748 


+ 




Hence the renormalization constant Zg of section |21 is given by 


9 

2^ ^ 2 


22 g^Ga M f f 


Z~^ = 1 

® ' 3e IOtt^ ' 3e IOtt^ ^ ' \^27e^ 81e J V IOtt^ 7 

and the corresponding /3-function is of the shape 
d 


fd{g) = ^ HZgg^) 


-1 


g 11 g^C 2 34 f g^C 2 Y 2857 f g^C 2 \ 
= —e - a -^- Q -:r- Q -^ 


3’^ 167r2 3 ^ V IOtt^ J 


54 ^ V 167r2 J 


(90) 


o{g^m 


agreement with former results 127103 Ei 


8 Conclusion and Discussion 

We have succeeded in computing the three-loop /3-function for a nonabelian 
gauge theory with a simple gauge group by means of the covariant background 
held method. Our results conhrm those of 03 and 03, where conventional 
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field theory methods were used. We also agree with the recent work of m 
who employed the noncovariant version of the background field method and 
extended the results of Em EE] to the case with Yukawa couplings and fermions. 
We are of the opinion that our work demonstrates that it is not pedantic or 
impractical to insist on manifest background field gauge invariance throughout. 
On the contrary, it leads to a considerable reduction in calculational labor. The 
extension to fermion and scalar fields is straightforward (see e. g. m for the 
exact propagator of a quark field) and would make our methods suitable for 
phenomenological applications. 

In the covariant background field method there are precisely as many vertices 
and propagators as in conventional field theory. We selected the Fock-Schwinger 
gauge for the background field because this field then appears only through its 
field strength. It follows that Feynman graphs are at most overall logarith¬ 
mically divergent. This leads to a considerable reduction in labor, especially 
in higher loop orders, as compared to the usual way of employing the back¬ 
ground field method where overall quadratically divergent graphs appear. An¬ 
other advantage of choosing the Fock-Schwinger gauge is, that in contradiction 
to [2niinni) the explicit knowledge of the Schwinger phasefactor implies factor¬ 
ization of Lorentz and Lie structure of the Feynman graphs. We showed that 
calculational labor can be further reduced by first using an identity relating the 
exact gluon and ghost propagators. Via such manipulations we could reduce the 
full two-loop calculation to a single logarithmically divergent Feynman graph. 
At the three-loop level, all tetrahedral graphs could be removed (see sect, d 
and hence simple poles with a residue proportional to C(3) do not occur in the 
first place. This improves on earlier observations of explicit cancellation of such 
poles m In fact, the efficiency of the covariant background field method allows 
us to use a high level language such as Mathematica to reduce the amount of 
Feynman graphs as described in section d 

We now give a brief discussion of the prospects for an eventual three-loop su¬ 
pergravity calculation along the lines of our work in Yang-Mills theory. We 
will use the component field language and comment afterwards on the use of 
superfields. It is well-known that supergravity gives rise to a one- and two-loop 
finite S-matrix m- However, at three loops one expects [SB] a divergence of 
the effective action of the schematic form 

= ^IC^ + ... , (92) 

where G is Newton’s constant, is a scalar consisting of four Weyl tensors 
and the dots represent fermionic terms. Only an explicit calculation can deter¬ 
mine the absence or presence of this potential divergence. Contrary to naive 
expectations, it was shown recently with string-theoretical methods m that 
N = 8 supergravity probably does not have such a three-loop divergence. In 
fact, the first divergence of maximal supergravity is expected to occur at five- 
loop order. It would be interesting to confirm this surprising result with an 
explicit three-loop calculation and investigate if it still holds for N < 8. We 
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shall now sketch how such a calculation would proceed by means of both the co¬ 
variant and noncovariant background field method. Starting from the classical 
action for supergravity, one would select suitable background covariant gauges 
for the various quantum fields (701 ini For the background field metric one 
should work in normal coordinates in [701133 [731, this being the analogue of 
the Fock-Schwinger gauge for Yang-Mills fields. Indeed, in this gauge one can 
express the background metric in terms of its field strength, i. e. the Weyl ten¬ 
sor. In addition we are free to demand a covariantly constant background Weyl 
tensor, i. e. we may assume a locally symmetric background space. 



Fig. a Fig. b 


Figure 12: Comparison between the ordinary (a) and the covariant (b) back¬ 
ground field method in gravity. 

In the noncovariant approach |28| one linearizes the background metric, = 
This ultimately leads one to compute three-loop graphs with maximal 
degree of divergence equal to eight, see figure IT^ . On the other hand, within 
the covariant approach described above, one must compute only logarithmically 
divergent three-loop graphs, see figure d) The exact propagators appearing 
in these graphs are of the generic form 

4 

G{x,x') = '^C'^Gk{x-x') , (93) 

fc =0 

where we must include terms of up to fourth order in the Weyl tensor in order 
to take account of all divergences. Supergravity has a non-polynomial action, 
so that vertices up to sixth order in the quantum fields must be taken into 
account, see figure 113 The graphs involving fifth and sixth degree vertices 
are however as easy to compute as two-loop and one-loop graphs, respectively. 
Thus, as in Yang-Mills theory, the most complicated graphs are tetrahedral. 
We therefore choose our gauges such that the number of three-quantum-field 
vertices is minimized. As was already shown in one can reduce the number 
of three-graviton vertices to just two. 

Of course, there also exist efficient superspace methods for such calculations, at 
least for TV = 1 or 2 supergravity. The background field method can be extended 
to superspace |73 and used for calculations at the one-loop level in supergrav¬ 
ity m and two-loop level in super Yang-Mills theories We should mention 
in particular m, which advocates the use of gauge-covariant superderivatives 
instead of ordinary superderivatives and employs exact super-propagators. This 
is shown to lead to improved power counting and hence less work. This is com¬ 
pletely analogous to our work in ordinary space. Just like the manipulations 
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Figure 13: Three-loop graphs in supergravity (thick lines represent all quantum 
fields to the extent that the corresponding quantum vertices exist). 


with gauge-covariant superderivatives directly on supergraphs, we could manip¬ 
ulate ordinary gauge-covariant derivatives on ordinary graphs. Selecting suit¬ 
able normal coordinate gauges in superspace as in m should lead to further 
improvement since only the background superfield strength would then appear. 


A Notation 

The greek letters ... are used to indicate lorentzian indices which can adopt 
values from 0 to d — 1 with d being the dimension of the coordinate space. 
The colored indices of SU(N) are indicated by the latin letters in the 

fundamental and by the latin letters a,b,... in the adjoint representation. 

Following DeWitt the dot product 

a-b = j dxa^{x)bi{x) (94) 

is defined in such a way, that i represents all indices, which are in common in a 
and b and have to be summed up. Furthermore we write dx instead of d‘^x, as 
long as the context allows no confusion. 

A bracket underneath a Feynman graph gives the value of the UV-divergences. 
They have to be read from the right to the left. With every comma the degree 
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of the pole part increases. The coefficients which are separated by semicolons 
belong to the same degree of divergences. They are only specified, if the diver¬ 
gences contain C(2n)-functions. The definition is given by 


(• ■ • ,o,5,ci;c2) 


a b ci-hC2C(3) 

Q “r o 


(95) 


B Identities of Gi and Ri 

In d dimensions the functions Gi are defined by 

G._i(x) = ^ for*>l (96) 

(47r)2 V 4 / 

Thus, for d < 2 all functions Gi are free of UV-divergences. This result can be 
analytically continued to higher dimensions. But in contrast for i > d/2 the 
functions Gi show up IR divergences, which means that as soon as in d = n — e 
dimensions the inequality i + \ > n/2 holds the functions Gi have to be replaced 
by the following functions 

— e / 2\ 

Ri-i{x) = Gj_i(a;)-f £- 2 ^ 44 j for i > d/2 , (97) 

where the added pole part cancels the IR divergence of the function Gi. 
Concerning the derivatives we get the following identities 


df^Gi{x) = 

-^x^Gi-i{x) 

df,Rt{x) = 

-^Xfj.Ri-i{x) 

d‘^Gi(x) = 

-i Gi-i{x) 

d‘^Ri{x) = 


d'^Goix) = 

-5\x). 


Using these identities, products of two functions Gi can be transformed in suited 
products. Some of the most useful transformations are listed below 


G^d^Gj 


G, d.d.Gi 


Gj-i df^Gi+i (98) 

(a^G,+i) a,G,_i - i <5^.G,G,_1 (99) 

(a^9,G,+2) G,_2 - i 5^, (G, G,_1 - G,+1 G,_2) , (100) 
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whereas no index adopts a value less than zero. Several identities depend on 
the dimension d. These are for example 


{d — 2 — 2i) Gi{x) df_iGj{x) = {d — 2 — 2j) Gj{x) df_iGi{x) (101) 

d^Goix)d>^G,+iix) = (^^-?jGoix)G,ix), ( 102 ) 

which have to be used carefully, since the value of d depends on whether it is 
part of an UV divergent subloop or not. 

The identities presented above deal with products of propagators depending all 
on the same endpoints which means that in coordinate space their values have 
just to be multiplied. 

Regarding propagators concatenated to each other we need the following prop¬ 
erty 

f , . 1 1 ^ ^d /2 nd/2 - a)r{d/2 - (3)r{a + (3 - d/2) 

J ^ {x - {y - r{a)T{(3)r{d-a-P) 


(a; _ z)'^{a+l3-d/2) ’ 


which leads to 


Gi{xo - Xi+i) = i\ j dxidx2-■ ■ dxiGo{xo - xi)Go{xi - X 2 ) ■ ■ ■ Go{xi - Xi+i) 


XQ XI X2 

= ^! a t t 


Thus the convolution of n -I- 1 functions Go multiplied by a factor n\ equals the 
function Gi. This explains why the graphical representation of the function Gi 
consists of a factor i\ and a line with i + 1 dots. 


C Group Theory 


For a simple compact Lie group G with elements 

U{x) = withA(a:) = taA“(x), 

where the antihermitian generators ta with a = 1,..., dim(G) satisfy the algebra 

— fabcicj 

the Dynkin index and the Casimir operator are defined by 

tr(t,4) = -T{D)Sab (105) 

tata = -G2{D)1 . (106) 



This representation is called fundamental and T{D) = 1/2 is a common nor¬ 
malization choice. 

Specializing on SU(N) groups the dimension equals iV^ — 1 and we get 


1 - 

ta ta = - 1 

2N 

(107) 

1 

ta tb ta — 2,N 

(108) 

Thus in the fundamental representation the Casimir operator is given by (1 — 
N'^)/{2N) multiplied with the unit matrix. Please note, that in Sab the indices 
run up to N'^ — 1, whereas in Sij the range of the indices is limited to N. 

Choosing the adjoint representation 


(A'a'jbc — fabc 

(109) 

we are lead to the following Casimir operator 


tviTaTb) = -NSab 

(110) 


for the SU(N) groups. 

In this paper all terms written in the adjoint representation are illustrated bold. 
The Casimir operator remains undetermined. The adjoint held strength e. g. is 
given by 


rpa (rria\bc jTia £abc 

= 

= Y'^pFp^ 

= TrF^^ 

= 

where tr respectively Tr are the shortcuts for the trace over the group indices 
respectively Lorentz indices. 

To calculate products of structure constants in an effective way we establish a 
graphical representation. Using the dehnitions 

5ab = a - b; Saa = dim(G') 

the table o shows some needed identities. 
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Table 3: Identities of the structure constants 


D Calculation of two two-loop graphs 

In this appendix we will calculate the UV-divergences of the last two graphs in 
the last row of figure El explicitly. Using 

+ F^aF,p5^5^G2 - ^F^^Gi + 0(f3) 
and 

- iF<,,F^,Gi - Var^p,dPd^G2 
+ F,,Fp^5^9,G2 + 0(f3) 
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we get 


Dfj,Ga-T D i, 



+ 2S,rF^xd^d,G, - 2d^rF,xd^d^G,) 

- - F,,Fp,dPd^G2 + F,rFppdPdM2 

+ d^rFpxF,pd^dPG2 - ^S.rF%Gl + ^S.rF^dpd, 

+ ^S^^Flpd^dPdpdMs - 2Fl^dpd,G2 
+ 0(f3). 


dpd,G2 


( 111 ) 


This result combined with the identity (l,r^i,,FpCT) = — (1,1, Fpi,FpCT)/2 and 
the symmetries of the 0 graphs leads us to 


+ 6 (Go, Go, GiF^) - (5p,Go, Go, J'^'GiF^) 
-^(<5m.Go,<5'^"Go,GiF2) 


-2((5p,Go,5'^"Gi,a'’a"GiF2^) 

-2(5p,Go,(5'^"9^a'^Go,G2F2J 


- 2 (Go, (5p.Go, 5^^’^dPd^G2Fl^) + 0(f3)(112) 


As expected, all graphs containing the field strength linearly vanish due to 
symmetry reasons. The first term in ins does not contribute either, since con¬ 
tracting the term d^duG^ with 5^" of the first propagator, we get a 5-function, 
which entails that initial and end point of the Go-function are identical. Since 
Go(0) equals zero by definition, the first term can be neglected. 

We refrain from contracting those 5-tensors, which lead to an additional factor 
of dimension d. Thus, when using the R* method, we can still decide, if all 
5-tensors building up a factor d belong to any subgraph and if therefor the e 
part of the factor d has to be neglected, while calculating the divergences of the 
appropriate subgraph. 

In many cases we can eliminate these 5-tensors with the help of the identity 


dG, = (2 + 2f-Xpa'^)G, 


(113) 


unless the Gi function belongs exactly to the same subgraphs as the 5-tensors 
do. Applying this identity e. g. to the second graph of (trnii . we get 
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whereas the last equal sign is based on the identity 

di-d2 = \dl-\dl-\dl (115) 

The numbers illustrate on which propagator of the 3-vertex the derivatives op¬ 
erate. In momentum space the analogue identity is known as 

PI-P2 = \pI - \pI - ^pI- (116) 

For the second graph of the last row in hgureElwe get, by using the identities 
and ifinnii . the following result 

= i(l,Gi,GiF2)+3 (Go,Go,GiF2) 
-4(Go,Go,5'^a'^G2F2J 
-2(Go,5'^a'^Go,G2F2J+0 (f3), (117) 

(Go, Go, GiF^) + (Go, Go, d^d''G2Fl,) 

+ 2(Go,9'^5"Go,G2F2^)+0(f3). (118) 

Next, we want to substitute the last two Feynman graphs of mg respectively 
both containing derivatives, with Feynman graphs containing no deriva¬ 
tives. For this reason we introduce two Feynman graphs without subdivergences. 
The first one is given by 


whereas the third graph yields 



{dpG^,d,G^,dPd'^d^^d''G2F%) 


+ 2 




i(l,Gi,GiF2) 

+ (Go,Go,9^9"G2F2J, (119) 



where we added the graphical notation to make things more comprehensive. It is 
easily shown, that this graph, by repetitive use of the identity (tTTCIl . decomposes 
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in two graphs. Both can be found in the equations itTTTI) and itTTHll . whereas the 
last one belongs to the group of graphs which we want to eliminate. Since the 
introduced graph contains no subdivergences, which means, that it has only a 
simple pole in e, we are allowed to use the identity (1^ with d = 4. Therefore, 
being interested in the divergent part only, the new graph simplifies to 



{dpG,,d,G,,dPd'^df^d‘'G2Fl,) 


-1 (1,1,F2)+O(e0) 

respectively 

-^{dpGud^GudPd^GiF^) + 0{e°) 

i(l,Gi,GiF2) 

-i(Go,Go,GiF2)+0(eO). (120) 


Finally, by comparing the right sides of the equations (I119II and (tnnii we get 
the identity 


{Go,Go,d>^d’^G2F%) = -i(l,Gi,GiF2)-i(Go,Go,GiF2)+0(e°), (121) 


which allows us to eliminate one of the said graphs. In a similar manner we 
create another graph containing no subdivergences 



= {dpGi,dPd^drG2,d'^d^d^^d''G2Fl,) ( 122 ) 


1/ 


which leads us to the identity 

(Go,9^9‘'Go,G 2F2J = -l(i,Gi,GiF2) + i(Go,Go,GiF2)+0(eO) . 

(123) 
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Taking in account the equation 


= 2(l,Gi,GiF2)+O(e0) (124) 


eventually only two graphs remain, which contribute to the divergence for the 
Yang-Mills-theory up to two loop as presented in sectional 


^2"" “ — {Gfiv^G paD u) — 2{G^inGpi,D u, ) 


, D„G„ 


^(1, Gi, GiF2) - ^(Go, Go, GiF 2) + 0(e°). 


(125) 


The result, which is wellknown from the literature, can be verified by looking 
up the divergences of these graphs given in appendix IeI 


E Divergences of some graphs 

In this appendix the divergences of some graphs are specified. Though, for the 
renormalization up to the order of two-loop, only the divergences of one two- 
loop graph namely 01.3 is needed, the other graphs are required while doing 
the three-loop renormalization, since they appear as subgraphs. 



00 01 



01.2 


02.1 02.2 



01.1 




01.3 


02.3 


01.4 



02.4 




Table El contains two results for every graph. The first column transcribed 
with K.R gives the divergences of the graphs with uncontracted derivatives, 
whereas the second column transcribed with ICRSpv gives the divergences of 
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Table 4: Divergences of some two-loop graphs 


the graphs with contracted derivatives. We want to emphasize once more that, 
since we work with renormalized graphs, the result of the second column cannot 
always be obtained from the first one by merely multiplying its result with a 
factor d. Again we have to watch carefully if the contracted derivatives belong 
to a subdivergence. 

Concerning the divergences of the graphs, table IeI shows, that the partial inte¬ 
gration holds, since for i = 1,2 the equations 

(0i.l)-b (0i.2)-1-(0b3) = 0 
2(0b2) -b (0b4) = 0 


are satisfied for both columns. 

The following identities do not apply to graphs with contracted derivatives. 
They are established with the help of some identities of appendix H and are 
given by 

( 01 . 2 )- 2 ( 02 . 2 ) = 0 
(4-d)(01.2)-(2-d)(01.3) = 0, 

whereas e. g. the equation 

(6-d)(02.2)-(2-(i)(02.3) = 0 

does not hold, due to the reasons mentioned above. 
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